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1 Introduction 

Let M be a compact smooth Riemannian manifold of dimension n with line 
element 

ds'^ = Qijdx^dxK (1) 

In order to investigate the Navier-Stokes equation on manifolds the Navier- 
Stokes equation system on flat manifolds 

f |^-j.Av+(vV)v = -Vp + fe., 

V • V = 0, (2) 
I v(0,.) = h, 

has to be reinterpreted. The class of flat compact manifolds where ([2]) is a 
correct formulation is rather limited (the classification is well-studied). The 
n-torus T" is the most natural example in the sense that this is the only 
type of fiat compact manifolds which occurs in any dimension n > 1. In a 
more general context of Riemannian manifolds it is natural to replace the 
nonlinear term by the covariant derivative of the vector field with respect to 
itself, i.e. 

(v • V)v is replaced naturally by VvV, (3) 



^ *This is an abbreviated version of a more extensive version to be submitted 
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where the latter symbol denotes the covariant derivative in standard invari- 
ant notation. That this requirement is consistent with the formulation on 
flat manifolds can be seen easily by writing the covariant derivative in co- 
ordinates. For the jth. component we have (we use Einstein notation in the 
following if convenient) 




k=l fc=l \ m=l 



where Cj denotes the j'th unit vector of the Euclidean basis. The Christoffel 
symbols become zero if the manifold is flat. Hence taking the covari- 
ant derivative is a quite natural extension which collapses to the classical 
equation term for the n-torus. Recall that 

= ^5^' {9jk,i + 9ik,j + Qij.k) , (5) 

where (g^') denotes the inverse of (gij). Another matter is the reinterpreta- 
tion of the Laplacian. We can deflne it in terms of the covariant derivative 
and its adjoint or in terms of the exterior derivative and its adjoint. The 
former possibility leads to the Bochner Laplacian 

Lb = -A*A (6) 

where 

A* : C°° (M, T* (g)T)^ C°° (M, T) (7) 

denotes the adjoint of the connection A on the tangent bundle on M and 
T = TM (resp. T* = T*M) denotes the tangent bundle (resp. cotangent 
bundle). This A is not to be confused with the usual Laplacian of course, 
but since we do not use this notation in the following there should be no 
peril of confusion). Another possibility is called the Hodge Laplacian which 
is defined by 

LH = -{d*d + dd*), (8) 

where d* denotes the adjoint. Note that Lh is related to a Dirac opera- 
tor D. Indeed from index theory we know that the difference of and 
the Bochner-Laplacian (applied to some vector field v) is given in terms of 
the curvature tensor of the connection on the tangent bundle (applied to 
v). The incompressibility condition usually leads to a simplification of the 
interpretation of the Laplacian, and in many cases it seems reasonable to 
interprete it by the Bochner-Laplacian plus the Ricci-tensor applied to the 
velocity field. In any case we may subsume a lot of possibilities by assum- 
ing that we may substitute the Laplacian by a linear scalar second order 
diffusion operator L which has local coordinates 

j,k=l ^ " k=l " 
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where the subscript U indicates that the operator is looked at with respect 
to a chart defined on [/ C M*^. WeU, the second order coefficients may 
also depend on a third index i, but let us keep things simple, since this 
is not essential, because we have no coupling of the second order terms in 
the Navier-Stokes equation. We may impose uniform ellipticity conditions 
on the second order terms, postponing possible generalizations involving a 
Hormander condition to subsequent investigations. All these local operators 
Lu together with an atlas on M define the global operator L on M which 
may be interpreted as a Hodge-Laplacian or Bochner-Laplacian in specific 
circumstances. In any case we assume that the operator in ([9]) is uniform 
elliptic with bounded smooth coefficients. Finally, we need to reinterpret 
the incompressibility condition on a manifold. Well, in Einstein notation 
(and with Einstein summation) this is just 

divv.= v^j. (10) 

The considerations so far lead us to the conclusion that we may define the 
Cauchy problem for the Navier-Stokes equation on manifolds by the equation 
system 

^ - z^Lv + Vvv = - VmP, 

divv = 0, (11) 
I v(0,.) = h, 

where V^./p denotes (some interpretation of) the gradient of the scalar func- 
tion p. On a manifold the gradient is naturally interpreted to be the eval- 
uation of a vector field at p. Locally, we can work this out as in Euclidean 
space, i.e. we may look at local coordinates in a certain chart with values 
in U where the Cristoffel symbols disappear and the coefficients gij of the 
line element satisfy gij = 5ij locally. We have to transform the coefficients 
Ojj and bi of the operator L accordingly. We call these coefficient functions 
'locally fiat on [/' and, keeping our general notation, we just write 

in order to indicate that we consider the coefficient functions in a chart 
where the metric is Euclidean and the affine connections involved collapse 
to a directional derivative. Then locally on U everything looks as in the Eu- 
clidean space (except that globally -as U varies -we have variable coefficients 
of first and second order, i.e., we have for all 1 < i < n 

9?^ _ v^n .J,U d'^v' v^w i.f,U dv' 
dt l^j,k=l"-jk dxjdxk l^k=l"k dxk 

(13) 

+Ei=iv:kv' = -p,i- 
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Note that the coefficients in (|12p are locally constant in locally flat coordi- 
nates in the usual interpretations of Laplacians on manifolds (however, we 
are flexible with regard to this point to some extent). This procedure is 
useful for elimination of the pressure from the equation. (We shall make use 
of nonzero Christoffel functions below.) Derivation for each i of equation 
(|13p with respect to Xi and summing up using incompressibility leads to a 
Laplacian equation on U in the form 

-Ap = 5(v,Vv)j,, (14) 

where 5-^'^ (v, Vv)^ : [/ ^ M is a function deflned by 

D {V, [X) .— 2^j=i 2^j,k=l "'jk,idxjdxk 1^1=1 l^k=l°k,i 

+ T.lk=iv]kV% 

(15) 

Again, the superscripts / and U indicate that we look at the functional S 
in a locally flat metric on a local domain U . Therefore for the first order 
spatial derivatives of the coefficients we have j(a;) = and 6^'- (x) = 
for usual interpretations of the Laplacian on the manifold M . In this case 
the equation (jlSp collapses to the usual Leray elimination of the pressure. 
This is especially true in the case of flat manifolds such as the n-torus. 

Next we have to discuss this local Poisson equation on [/ C W^. Actually, 
in terms of a given vector fleld v we have a Dirichlet problem here of the 
form 

Ap(t, .) = / on [/, 

(16) 

p(t, x) = g on dU, 

where dU is the boundary of U. Note that we have a family of Dirichlet 
problems here (one for each time t > 0), where t serves as an external 
parameter. In our case / and g are defined in terms of the same function 
S (v, Vv) , where for fixed t and x & U 

/(x):=-5^'^(v,Vv)(t,x), (17) 

and 

5(x):=-S/'^^(v,Vv)(t,x). (18) 

Here, we may assume that domain U is small enough such that the definition 
of S^'^ in flat coordinates can be extended to the boundary of dU of U. 
Another possibility is to deflne S^'^^ (v, Vv) (t, x) with respect to a finite 
set of charts with image V^, i G / such that dU C Ui^jVi such that for fixed 
t and X G dU n we have 

gf,aunv, (v, Vv) (t,x) = 5^' (v, Vv)^, (t,x), (19) 
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where the charts are identical on U OVi for all i £ I. Here and in general if 
we consider the functional S" in a certain unspecified chart on [/ C M"" then 
we just write to indicate the functional S in local coordinates at U. We 
may assume that the atlas for the smooth compact manifold M is chosen 
such that each Ui has a C°°-boundary S. As it is well known the solution for 
the pressure p on the domain U can be obtain by the sum of two solutions 
pi and p2 of two Dirichlet problems with zero boundary condition for pi and 
zero interior condition for p2 such that 

P = Pi+P2 (20) 

where both functions are given in terms of the Green's function Gjj (here 
again the superscript / indicates with respect to locally flat coordinates), 
i.e., 

Pi= [ f{y)G(j{x,y)dy, (21) 
Ju 

and 

P2= [ g{y)d,yG{j{x,y)da{y), (22) 

JdU 

where dU denotes the boundary of U and d^^ denotes the derivative with 
respect to the outward normal. Hence locally on U and with respect to 
locally flat coordinates we have a Leray representation of the Navier-Stokes 
equation of the form 

/[/ m-G{j{^^ y)S^''' (v, Vv) (t, y)dy + /^^ d,yG{j{x, y)Sf^^ (v, Vv) (i, y)dy. 

(23) 

In general local coordinates we get the same equation in the form 

dt 2^j,k=l '^jk\-^) dxjdxk l^k=l "k y-^l dxk 

+ ELi {v]k + Tl=iv'^'^'fk) = 

lu £rGu{x, y)S^ (v, Vv) {t, y)dy + j^^ d^^Guix, y)S^ (v, Vv) (t, y)dy, 

(24) 

where the superscript U to the coefficients aij, hi, and F*^ indicate that this is 
the local representation in a unspecified chart -0 : Um C M ^ U. Note that 
in (fM|) we dropped the superscript / (indicating locally fiat coordinates). 
The related terms are considered to be defined in terms of the terms with 
superscript / via coordinate transformation (this defines af, S'^ etc.). Now 
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if we take a finite atlas of the smooth compact manifold M of the form 
A := {ipm ■ Umt^ ~^ ^m, m € 1} (where we assume that each Um has a 
C°°-boundary such that the Green's function Gjj^ exists, then we obtain a 
representation of the Navier-Stokes equation in Leray projection form on M 
by a family of local equations each of which is of the form ()24p . Note that 
except for the term 

/ ^,^Gu{x,y)S{^^,\/^^){t,y)dy (25) 
JdU 

all terms in (|24|) are local. Indeed the term (|25|) describes the communication 
for a system of local equations that we are going to define. 

Next we derive this global system of local equations of the form (|24p using 
a partition of unity. Using this partition of unity we can derive localized 
equations where the fiow of information over the boundary (the coupling of 
the local equations) is based on the terms (|25p . First, let us recall some 
facts concerning partitions of unity on smooth manifolds. We have 

Proposition 1.1. Let M be a compact manifold and let (Wi)jg/ be an open 
cover of M. Then there exists an open cover {Vj)ji^j subordinate to {Wi)i^i 
and a family of real smooth functions such that 

i) (pj : M , where supp{(j)j) C Vj, and 

ii) 4>j{x) > for all x E Vj, and for all x ^ M 

Y.H^) = l (26) 

Remark 1.2. Since M is compact we may assume that the cover (V^Ojej is 
finite, i.e., we may assume that the index set J is of finite cardinality. 

Remark 1.3. The above result may be generalized to manifolds which are 
countably unions of compacts. In this paper we restrict our considerations to 
compact manifolds. This simplifies some steps of the proof of global existence 
which are more complicated in the context of non-compact manifolds. 

We apply the latter proposition and consider a partition of unity ((/)j)j£j 
subordinate to (yj)jgj where we may assume that the latter cover is such 
that we have a description of the Navier-Stokes equation in locally flat co- 
ordinates on each V^-, i.e., for each Vj we have a chart ijjj : Vj — > tpjiVj) =: 
Uj C M", and such that the equation takes the form p^ . or, equivalently, 
the form ()24p in the chosen coordinates. We shall define a system for the 
family (i'*"');^<j<„ jgj) where f*-' denotes the restriction of to the domain 
Uj in a chart ipj. Note that for all j € J we have in local coordinates (on 
Uj) 

v^^{t,x)=^<P%ixy{t,x) (27) 

k&Jj 
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along with 

Jj:={keJ\UjnUk^0}, (28) 
and where for ah j,k J we define 



ct>k [ipjHz)) if i;f{z) e Vk (29) 
else. 



The solution of the system for the family of functions (^*"')i<j<„ gj provides 
enough information for reconstruction of the solution v of tTie Navier-Stokes 
equation on a manifold. The family (f^*"''') i<j<,„ jgj along with 

v'^^ : [0, oo) xUj 

(30) 

The latter family, as derived from the family ■yi<j<„ jgj may be considered 
as a representation in local coordinates of the solution of the incompressible 
Navier-Stokes equation on a compact Riemannian manifold subordinated 
to a finite atlas which allows for locally flat representations. In locally flat 
coordinates we get an equation 



(31) 



where 



+ ELi = ^Lt. + ^coup. (v, Vv) , 

^int,i = X ^^^^ ^^^^ 



along with 

(V, Vv) (t, X) := - EILi .=1 «^L(^)^ 



dXmdXk 



(33) 



■spn v^n ('r\ "'" I 



and the coupling term "S^Q^p ■ (v, Vv) is defined below. Recall that in a 

usual interpretation of the Laplacian on manifolds in locally flat coordinates 

f u 

the spatial derivative of the second order and first order coefficients aj^ ^ 

and fe^'- are zero and the definition in (|33|) simplifies accordingly. However, 
tracking them has the convenience that we can easily retrieve the information 
how the solution may look like in general coordinates. Moreover, we cannot 
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stick to the same locally flat coordinates for v^^ if we consider the functions 
u*P on the adjacent domains Up with p £ Jj\ {j} in general. It is natural to 
consider them in natural local coordinates, i.e. such that the line elements 
looks locally like ds^ = "^^p gmpdxmdxp. Then we have to substitute usual 
derivatives by affine connections. However, on the boundary of dUj we may 
consider locally flat coordinates (w.l.o.g. we may assume that the partition 
of unity is fine enough to have flat coordinates in a small neighborhood of 
each coordinate patch Uj). Next we observe that 

v'^i = 4>^-v'^ (34) 

has zero values on the boundary dUj . Accordingly, in order to define S^Q^p ^ 
(which denote coupling terms) we consider for p G Jj\ {j} and x £ UpH dUj 

Z^i=l Z^fe=l "fc,j y-^J dxk "I" Z^i,fc=l ^fe 



(35) 



where the notation a^^'''^^^^ and l)^^^^^^^^ indicates that these are coeffi- 
cients which are flat in C/j-coordinates (but not necessarily flat in coordinates 
which are flat in Up). Hence, in a usual interpretation of the Laplacian on a 
manifold in locally flat coordinates theay are locally constant and their par- 
tial derivatives are zero. However, denoting these partial derivatives in ()35p 
accounts for a more general situation (possible generalisation) and reminds 
us how the term may be rewritten in general coordinates. Next we have 

•5coup,i (v, Vv) = 

(36) 

Epgj,\0} hu d.yGuix, y)5/'^P'^ (v, Vv) (t, y)dy. 

Again, we emphasize that the boundary terms for j = p disappears since 
(j^'jj^-^ becomes zero on the boundary of Uj. Now we have a system of card( J) 
equations of form ([HI]) where equation number j is coupled with equation 
number k ii k £ Jj . 

In general coordinates we write Sj^'^ = S^'^^'^ , i.e., we drop the super- 
script /. In this representation the derivatives if the coefficient functions 
across the boundary dUj are no longer zero in general. Indeed, they de- 
scribe some part of the flow of informations or the coupling between the 
local Navier Stokes equations for v'^K 

This leads to the idea that we may define a time-local scheme by solving 
equations for v"^^^ which approximate v^^ for 1 < i < n and j € J, and 
where we take the coupling information from the previous iteration step 
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m — 1. More precisely, for j(zJ,l<i<n, and m > 1 define recursively 

at l^j,k=l "■jk \-^) dxjdxk 2^k=l "k \-^) dxk 

+ EUvrv'^- (37) 
= 5^;;;;jv,Vv) + 5eX"'(v,Vv), 



where 

^inT. = / i^.^U. (v, Vv) (i, y)dy, 

JUj i 

along with 

sr"" (V, Vv) (t, X) := - Er=i Yi,u=. 

1^1=1 l^k=l"k,i y-^) dxk ~^ ^i,k=l",k",i ' 

and 

•^coup'^ (v, Vv) = 

SpeJAO} iaf/ ^-.G't;(rE, y)5/'ii'.™-i'^ (v, Vv) (t, 2/)dy, 

along with 

SS.V^-^^ (V, Vv) (t, .) := - EILi E^,.=i 

_v^n. v^n ^/.'^p/^ m9i>'J("'-^) , y^n. ij{m-l) kj{m-l) 

2^1=1 2^fc=l V-^^ dxk 2^i,fc=l '^,fc 



(38) 



(39) 



(40) 



(41) 



We may start this local scheme at some time point t = to with v(to, •) = 
{v^{tQ, .),••• , v"'{tQ, .)) , if this function is known. Furthermore, we introduce 
boundary conditions in order to make this time-local scheme spatially global, 
i.e. for some ti > to, and for all j G J we add the boundary condition 

v'^"'it,x) = ^^^■^•(™-i)(t,x) for ah {t,x) G [to,ti] x dUj. (42) 



For m = 1 we may set v'^^^{t,x) = v^^^{t,x) := v'''{to,x) for all {t,x) G 
[to,ii] X dUj. Especially, if t = and m = 1 , then we define for 1 < i < n 
and j £ J 

^ij{m~i)|^^^ = v'^^ = h'^ = h'luy (43) 

For a horizon ti — to small enough we shall see that this scheme leads to a 
time-local fixed point iteration in classical space which defines in the limit 
m t oo a spatially global amd time-local solution 

V :[to,ti]x M ^TM (44) 
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of the incompressible Navier-Stokes equation on the manifold M. In a first 
step of our proof of global regular solutions we shall show that a refor- 
mulation of the scheme defined by the equations (f37|) . and (fl3j) and with 
boundary conditions (j42p converges to the local solution of the Navier tokes 
equation which is valid for small time t > 0. Note that the local right 
side of the Poisson equation which eliminates the pressure involves terms of 
the from ^^jV'^jiP^ (in locally flat coordinates), where the difference of two 
consecutive iteration steps (indexed by m) involves localisations of 

s-^n „,i(m-l) j(m-l) sr^n i{m-2) j{m~2) 

(45) 

along with 6v'^^^ := w"'/"* We shall use this in order to prove 

local convergence. 

The incompressible Navier-Stokes equation cannot be solved by a simple 
global fixed point iteration. Similar as in the case of the multidimensional 
Burgers equation and in the case of our global scheme for the incompress- 
ible Navier-Stokes equation on the whole domain of M" we choose a time- 
discretized scheme and construct fixed points which are local in time (cf. 
[1] and references therein). First note that it is convenient to choose a time 
step size < p < 1 and solve at each time step / > 1 for 

^P,i ^ . . . ^yP,i,n^^ : [/ - 1,/] X Af ^ TM (46) 

the incompressible Navier-tokes equation on the domain [/ — 1,/] x M in 
transformed time coordinates r with 

t = pr (47) 

and with initial data v''''(Z — 1, .) = vf'^~^{l — 1, .) being the final data of the 
previous time step number / — 1, where v'''^(0, .) = h(.). Here, as usual, the 
symbol TM denotes the tangential bundle of the manifold M. 

Remark 1.4. Note that in (j47p we choose a time step size p > which is 
independent of the time step number I. This is a difference to our scheme 
in [1] where we chose a 

Pir^J (48) 

The method of control introduced in this paper may also be used in order 
to simplify the argument on infinite domains. 

However, a more interesting fact about our scheme of a Navier-Stokes 
equation on compact manifolds is that we can operate on classical function 
spaces. We shall show that we have classical solutions for the family 

v'^ = v'luj, 1 <i <n, j e J, (49) 
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where the related family 

v'^^ = ^jjv'\ l<i<n, j eJ (50) 

may be considered as a representation of the solution of the Navier Stokes 
Cauchy problem on compact manifolds. Note that (09]) can be reconstructed 
from (|5Up and vice versa. Next we introduce a classical function space 

Proposition 1.5. For open and bounded C M" and consider the function 
space 

(Q) := | / : J7 ^ M| 9°/ exists for \a\ < m 

(51) 

and d^f has an continuous extension to J7| 

where a = (ai, • • • ,a„,) denotes a multiindex and 9° denote partial deriva- 
tives with respect to this multiindex. Then the function space C™ (fi) with 
the norm 

\f\m := |/lc^(n) := E \^"f\ (52) 

\a\<m 

is a Banach space. Here, 

\f\ :=sup|/(x)|. (53) 

Local convergence of the scheme for t;*-?"^ above implies that for given 
t > and for each j G J a, sequence 

yijM(^t^ .) = (j)UyVrn^^^ g C^{Uj) (54) 

converges in C'^{Uj). This leads to classical solutions v'^^^ in local spaces 
C^'^ {{tQ,ti) X Uj). We shall construct a bounded classical solution 

V : [0, oo) X M TM, 

where all components Vi are globally bounded with respect to the |.|i^2 norm, 
i.e., the suprema up to first order time derivatives and up to second order 
spatial derivatives are bounded in all local charts of a smooth atlas. 

In order to make the scheme global we introduce another idea, a control 
function r. First we introduce a new time variable 

r pt, (55) 

where /? > is a small constant to be determined, and such that each time 
step with respect to the new time-variable r is of unit length. Accordingly, 
for each time step number I > 1 the original velocity function restricted to 
the time interval [1 — 1,1) and denoted by v' = {v^'^, • • • , t;''")^ along with 

: [p{l -l),pl)xM ^ TM, (56) 
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there exists a time-transformed velocity 

vP'^ -.[1-1,1) X M ^TM (57) 

with components ti'''''*, 1 < i < n. Next, for integers I > 1 (time-steps) we 
shall construct a family of recursively defined functions r', Z€N, 1 < i < n, 
where 

-.[1-1,1] X M ^TM, (58) 
and consider for 1 < i < n and j J the equation system for the functions 

The series (r^)^ defines a global control function r which looks locally on Uj 

as a n-tuple (r''-^-', • • • , r^'"-')^ on each subdomain [1 — 1,1] x Uj (in a chart). 
These local representations determine the global control function 

r : [0, oo) X M ^ TM, (60) 

which is designed in order to control the absolute value and the first deriva- 
tives of the functions \r^'P'^ and of the functions r' themselves. We use the 
freedom we have in order to define the control function r''*-' in order to 
determine source terms of the equations for v'''^'', where we define r' by 
dynamic recursion such that r is a bounded function. If we know that r is 
bounded and globally Holder-continuous and has bounded continuous spatial 
derivatives of first order, and if we know that = v + r is globally Holder 
contiuous with bounded continuous derivatives of first order, then v = v*" — r 
is globally Holder-continuous with bounded derivatives of first order. Then 
we may look at the original Navier-Stokes equation and consider the first 
order coefficients involving the velocity components f * of n equations to be 
known coefficients of n linear scalar parabolic equation with a source term 
involving first spatial derivatives of i)* which we can consider to be known, 
too. Classical representations of the velocity components in terms of a fun- 
damental solution r„ are then available, and this leads to classical regularity 
of the velocity components v^. The idea for the construction of the control 
function r is as follows. We define the function r' for / > 1 time step by 
time step where at each time step certain source terms (consumption terms) 
are defined in terms of the data — 1, .) and r'~^(/ — 1, .) obtained 

at the previous time step. Now if we write down the equation for v'' we get 
bundle of Navier-Stokes-type terms for v'' plus a bundle of Navier-Stokes 
type terms for r plus mixed terms which are bilinear in v*" and r. Given 
v^''~^(Z — 1, .) and r''~^{l — 1, .) at time r = / — 1 we may determine r- via 
linearized equations with first order coefficients r^~^, and with some con- 
sumption terms or source terms. Next we define this controlled scheme in 
detail. 
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2 Definition of the controlled global scheme 



As indicated in the introduction, we shall use the time coordinates t = pt 
for some small p to be determined. This leads to a factor p for all terms 
except the time derivative if we replace t by r. At each time step / > 1 we 
have to solve for n ■ card (J) equations for 

yP^^^'^ : [1-1,1] xUj ^R, (61) 
for 1 < i < n and j € J, and where 

vP'^''^{t,.) =v^''^{t,.). (62) 
This family of local functions determines a spatially global function 

vP^^ :[l-l,l]x M ^TM (63) 

at each time step I > 1. The superscript p indicates that we are considering 
time coordinates r related to a time-step size p, and the number / > 1 
indicates the time step number. The equation for yP'^'^^ is based on ()3ip and 
is of the form 



(64) 



At each time step I > 1 the initial values of the functions v^'^^^{l — 1, .) are 
the final values of the previous time step, i.e., 

yP,i,ij{l - 1, .) = vP'^-^''^{l - 1, .). (65) 

Note that the local scheme in (|37|) represents a family of local Navier-Stokes 
equations in Leray projection form with an additional coupling term. In 
order to solve this system for each m > 1 we need an additional itera- 
tion. At each iteration step m we first the coupling term using the results 
yij{m-i)^ I < i < n, j G J from the preceding iteration step (m — 1) and 
then we solve iteratively linear equations at each substep p > 1 for functions 
y^jiTT-p approximating i;*-?™. The latter solve equations of the form 

dr ~ P 2^q,k=l^gk \^ ) dxqdxk ~ P l^k=l'^k dxk 

+P ELi vP^^''^V'''^"'^P-^^ = pS(^^ . (^yP,idMp-i) ^y^p,i,jHp-i)') (66) 
+P^coup,i (v/''''*^^™-!), Vv/''''*J('"-i)) , 
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where for p = 1 we have t;P.'.'=i"^(p-i) = yP^hkjm{o) ._ ^p,i,kj{m~i) _ 
following we also write 

^iiu"'^^ (V, Vv) := 5^^^ ^ (^^P^l^Jrnip^^) ^^^p,l,jmip-l)^ (g^^ 

and 

if this is convenient. Note that we have linearized the convection term and 
the 'trivialized' the Leray projection term (v, Vv) in the sense it 

is defined in by the previous iteration step {p — 1) and serves as a source 
function of a linear parabolic equation. At each (sub-)iteration step p we 
define 

^P,i,ijmp^l - 1, .) = vf'^-^''^l - 1, .), (69) 

of course. Furthermore, for each j G J we add a boundary condition for 
the local problem on [/ — 1, x Uj such that for all (r, x) € [Z — 1, /] x dUj 
the restriction x9C/, of w^-'.^J^p to the boundary [I -l,l]x dUj 

satisfies 

^^•'•^^■'"^l[z-i,z]xai/,(r,x) = vP^'^'^^^^-^\t,x). (70) 
Here recall that 

yP,l,ikk{m-l) _ ^U^p,l,ik{m-1) ^ j^y-j^-j 

and that Jj is defined in (j28|) Note that we use the partition of unity here 
in order to ensure that the latter prescription ()70p is well-defined for each 
iteration step m > 1. We shall see later that iteration with respect to p and 
then with respect to m leads to the time-local condition that in the limit for 
all j, k,q ^ J and all (r, x) € Uk Ci Uq D dUj ^ we have 

vP'^''''{t,x) = vf'^''''{T,x). (72) 

Furthermore, in case m = 1 we have to supplement 

v'^'^''^"'-'%-i,i]xau,ir,x) = vP^'-'^^\t,x). (73) 

for all (r,x) £[l-l,l]x dUj. 

The first step for a global existence proof of classical solutions of the 
incompressible Navier-Stokes equation on compact manifolds is to show the 
time-local convergence of this scheme. 

The next step is to define a global scheme. First we shall write down 
the system for v^'^'' := v^'' + r' on [I — 1,^] x M. Then we shall define r' 
via a linear equation with a right side (f)'' which involves the data v^'^''~^(Z — 
1, .) and r'~^(/ — 1, .) obtained at the previous step. Later in the proof of 
convergence we plug in our equation for r' into the right side of the equation 
for v^'P''- in order to show that the growth of v^''''' is uniformly bounded 
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independently of the time-step number / > 1. The construction of the 
iteration scheme describes the road on which we proceed in order to prove 
the global existence. We shall comment on this later on. First the equation 
for :[l -1,1] xUj follows from §6^ via the definition 

yr,p,l,ijmp _ yP,l,ijmp _|_ ^ ^j^-^ 

Indeed the family of control functions r''*-', 1 < i < n, j € J will be 
determined first by a linear equation in a first substep of time step / and is 
fixed then. Therefore, r''*-^ bears no iteration index m and no subiteration 
index p. We shall define the equations for r''*-', 1 < i < n, j € J below. 
The natural iteration scheme for is 

- dr P l^q,k=l'^qk dxqdxu P2^k=l"k \^)dxk I'^J 

+P ELi rfr^''^ + pSli^;r~' (V, Vv) + pS^i'^^-J"^ (v, Vv) , 
where 

Sr^P'-' (v, Vv) = S^^-' (v + r, Vv + Vr) , (76) 
^coupT"' (v, Vv) = 5^0'-;'-^ (v + r, Vv + Vr) (77) 

and where 

yr,p,^jmp^l - 1, .) = w'^'/'''-^'^J'(/ - 1, .). (78) 

Again, for each j € J we have boundary conditions, i.e., for all (r, x) S 
[l-l, I] X dUj the restriction x9c/, of to the boundary 

[Z — 1, X dUj satisfies 



^^''''''^^■"^^|[z-i,z]xat/,(r,x) = v'^''^'^'''^"'-'\r,x). (79) 

k&Jj 



In case m = 1 we define 

v''^''^'^'^^"'-'%-i,i]^dU,ir,x) = v^'P''-'''\t,x) (80) 

for all (t, x) G [/ — 1,/] x dUj. The global scheme is an iteration in time of 
a local iteration scheme of this controlled equation and a definition of the 
control function r'. At each time step I > 1 the functions v^'^'^~^'^^ {I — 1, .) 
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and r' ^'*-' (/ — 1, .) are known for 1 < i < n and j S J. For Z = 1 we may 
set r'"^'*-' = and 

^r,p,0,ij(/ _!,.)= i;P'O.^J(/ - 1, .) = h'^{.). (81) 

Next we define r''*-'' by the following n ■ card(J) equations. Note that the 
terms u^'/'''"!'*-? and r'"^'*-' in (|82|) below are abbreviations for the same 
functions evaluated at Z — 1, i.e. they equal v^'''''^~^''^^ {I — 1, .) and r^~^'^^^{l — 
1, .). The defining equation for r''*-'"^ (iteration index ni is in order to define 
a spatially global solution of this linear equation on the whole manifold M) 
is: 



4-/1 J'^i'^J-l,kjm 
'^Pl^k=l',k ' 

where we have to define Here, the upper index I — 1 indicates that 

we evaluate related functions at data obtained from the previous time step. 
Especially, 

^[nt'i '^"^ = "^int i (v^'^'^-^'^^"", Vv^'^'^-i'^J"^) . (83) 



(82) 



The initial conditions for (|82|) are 

r''^^"^(Z-l,.) =r'-i'^J(Z-l,.). (84) 

Since we want to construct a global solution 

r' : [/ - 1,/] X M ^ TM (85) 

to this linear parabolic equation, we have to ensure that the local solutions 
match on the boundaries dUj , j G J. We ensure this by boundary conditions 
in the same spirit as before, i.e., we define for all t,x) & [I — 1,1] x dUj 

r'.*i'"(r,x) = r'''*^''^'""^^(r,x). (86) 
keJj 

We shall define in terms of v^''P''-~^'^^ {l — l, .) and r'~-^'*-'(Z — 1, .). Note 
that 4>^'^^ is the only term on the right side of the equation ([82]) without a 
factor p (the time-step size) . For this reason these source terms can dominate 
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all other source terms which are approximations for corresponding terms in 
the equation for v^'^''"'^^. We define in terms of a sum = c^'''''*-? + 
. For each 1 < i < n and all j €z J the first summand c^^'''*-' is designed 



x,p,l,ij 



m regions 



in order to control the growth of r ''-^ and the growth of v' 
where A^~^^^^ (l — l, .) or first spatial derivatives of this function are large, and 
the summand cj)^''''^^ is designed in order to control the growth of v^'^''"'^^ in 
regions where v^''''''~^'^^ {I — 1, ■) or first spatial derivatives of this function are 
large and r^^~^'>'^^ [l — 1, .) is small (and therefore cannot control the growth 
of the function . There are several possibility to realize this. However, 

it seems to be most simple to put the growth control via — 1, .) on a 

different scale, i.e., to have an additional factor ^ in the definition of (f>^'^'^^ 
along with C S> 1 compared to the definition of cff'^'^K This way the growth 
control of v^'^'^'^^ does not interfer with the growth control of r''*-' in regions 
where the latter function becomes large. Note that the small time-step size 
p > makes it possible to control the growth in a region of space. The 
diffusion effect in the iteration of scalar parabolic equations is global, but it 
is small outside certain local regions because of the small time step size. We 
can control this by Gaussian estimates. We shall choose C ^ 1 such that 



J^ii I 



< C, and \v 



r,p,l,ij I 



and for all 1 < m < n 



Then we may define 



and 







Qyr,p,Lij 


< C, and 


dXm 





< c 



< c. 



(87) 



2_^r,p,(l-l),ij 



(x) 



(89) 



(90) 



^ '(x) :— ^2 

Note that we put ()90p on a different scale. 

The constant C ^ 1 will be chosen below. Note that we may choose a 
atep size of order 

(91) 

This means for 1 < i < n and j ^ J that both functions 0^'''*-? cf)^''-''^^ can 
dominate all terms which have a factor p in critical areas of growth. Next to 
dimension it depends only on the initial data h, constants of the manifold 
(determined by the Christoffel symbols or by the Riemann tensor) and the 
viscosity v > 0. Now that we have defined the source function for 
1 < i < n and j € J we can formulate the complete scheme. At the first 
time step / = 1 we set r^'*-' = for all 1 < i < n and j S J. Note that in 
this case 

(92) 



V 



r,p,0,ij 
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for all 1 < i < n and j € J. At each time step we define a spatially global 
iteration with iteration index m which has a subiteration with iteration 
index p. For each m > 1 we solve card(J) local Navier-Stokes equation via 
iteration in p, where 

dr Pl^j,k=l"'jk y-^J dxjdxk P l^k=l"k \-^) dxk 

+0 y^^PXiimp r,p,l,kjm{p-l) 

(93) 

with initial data 

^r,p,Mimp^Q^ _) ^^ijg^ (94) 

At this first time step for the subiteration with iteration index p we have 
the boundary condition 

v'''''''''"'%,i]>cdU,{r,x) = J2 i)(r,x), (95) 

keJj 

where the double superscript kk on the right side of (jl80p indicates the 
involvement of the partition of unity as explained above. 
For a step size p > which is small enough the sequences 

(^.,P,i..-P)^^^ (96) 

converge classically to a limit 

for all 1 < i < n and j S J . For the functional sequence (''^'^''''^'*"'"^)^gpj we 
get another spatially global iteration scheme 



dr Pl^j,k=l"'jk y.-^) dxjdxk P l^k = l"k y-^l dx 

n r,p,l,ij7n r,p,l,kjm 



V 



(98) 

with initial data 

v'''P'^^'^"'{0,.) = h'^{.), (99) 

and boundary condition 

v'''''''nm.dU,ir,x) = v^^p^'^'^'^"^-'\t,x), (100) 
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and where for m = 1 we have v'^'P'^'^^^'^ ^) = v''''''^''-'^'^) := h. 
For a step size p > which is smah enough the sequences 

K''''"^'")™.N (101) 

converge to a classical limit 

ptoo 

for all 1 < i < n and j G J. We shall see that for small p > this is a time- 
local fixed point iteration which leads to a spatially global and time-local 
solution 

v''''''^ : [0, l]xM TM. (103) 

This describes the first time step. Recursively, if v'^'P'^~^'^^ , and r'"^'*-' are 
known, then we construct (j}^^ as decribed in ()89p and ()90p above. Then we 
determine r''*-' via the linear Cauchy problem (|82|) . All this is done for all 
1 < i < n and all j G J, of course. The iteration for the control function 
j,i,ijm j^g^g ^ global limit which is locally r''*-^. The spatially local iteration 



scheme for the sequence (^v'^'P'^'^^^^p^ 



p6N 



IS 



Q^r,p,l,ijrnp flJ. ^ q2 ^r,p,l ,ijm.p uf'^t \ 9£[:^^^':^£™f, 

Or ~ P 2^j,k=l"-jk dxjdxk ~P2^k=l"k dxk 

1 V^jf''''^^'^^V^'P''-'''^^'^P~^^ 

= 4>^''' -pY1=i rfv-'^P^^-^'^^ -pEk=i vj'^-^''^r^-^'^^ 

-P ELi rfr^-^'"^ + p ELi rfv^'f'^'^^"^^P-^) + p ELi ^ J'''^^^^^"' V'''^^' 

+PEUrfr^''' 

-L^C-? (^r,p,l,ijm{p-l) y^r,p,l,ijm{p-l)\ 

^ int,i ^ ' ^ 



(104) 



+/0'5coup,i (v'''''''*J(™-i, Vv'''''''*J(™-i)) , 
with initial data 

^r,p,^jmp^l - 1, .) = v'^'P'^-^'^^il - 1, .), (105) 

and the boundary condition 

v'^'^'^''"'%,i]xau,ir,x) = ^;^'^'''*''('"-')(r,x), (106) 
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Here we plugged in a relation defined by the equation for the global control 
function r''*-' (the limit of r''*-''" as m t oo). For a step size p > which is 
small enough the sequences 

r,p,l,ijmp\ (107) 



converge classically to a limit 

p\oo 

for all 1 < i < n and j € J. 

Finally, for the functional sequence iy^''^'^'^^^') ^iz^ g^t another spa- 
tially global iteration scheme 

g^r,p,l.ijm ^ n f,U/ \ Q2^r,p,l,ijm ^ hf^U/ \ g.ur,p,l,ijm 



k 



r 



, v-^n r,p,l,ij(m—l) ; fc,- , v-^n l.ij I ki 

"^•^int'i ~ ^-^ciupj (v, Vv) 

(109) 

with initial data 



V 



r,p,l,ijm 



{l-l,.)=v''^P^^-^^'^{l-l,.), (110) 
and boundary condition 

v'-^'^'^'n[o,i].dU,ir,x) = J2 v''''''^""'^"'-'\r,x), (111) 

We shall see that for small p > this is a time-local fixed point iteration 
which leads to a spatially global and time-local solution 

v*"'^'' -.[l -1,1] X M ^ TM. (112) 

The following argument is essentially constructive up to the point that 
the choice of the time step size p > and the constant C > may be 
analysed more constructively. This will be done in a subsequent paper. At 
this analytic stage it suffices to sow that the constants p and C depend only 
on the data h the viscosity z^, and the coeffcients Qij of the line element of 
the underlying manifold M. 
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3 Main theorem 



We shall assume 



h G C°° (M, TM) 



(113) 



or, in a family of local charts ■0j '■ ~^ Uj C M", j (z J covering the 
manifold M we have components h^^^ € C°° (Uj) for all 1 < i < n and 
J € J. In this paper we set external forces to zero, although there is no 
problem to include them into the scheme we propose. This is just for the 
sake of formal simplicity of the description. An equivalent assumption on 
the initial data h is that the functions hi are located in Sobolev spaces of 
arbitrary order s G M, i.e., for all s € M we have 



This is the assumption which we used in the case of an infinite domain (cf. 
[T] ) . As indicated in the introduction our proof of a bounded regular solution 
of the incompressible Navier-Stokes equation consists of three main ideas: a) 
we introduce a time discretization and a series of linear time transformations 
t = pT such that time step size 1 in r-coordinates is related to a small time 
step size in original coordinates and small coefficients of spatial derivatives 
in transformed time coordinates such that an iteration procedure leads to 
a local solution in time. The choice of a constant step size p > 0, i.e., a 
step size which is independent of the time step number I is a feature which 
is possible in the context of compact manifolds M; h) in order to control 
the growth of the solution we introduce time-step by time step a control 
function r and solve an equivalent problem for the function v^' := v + r at 
each time step where r is itself a bounded function. The function r solves 
a linearized Navier-Stokes type equation, but with a 'consumption' source 
term which has been explained in the introduction to some extent. At 
each time step this source term consists of n components = ■ ■ ■ , 4>l^) 
such that the functions v''''''' and r' are bounded independent of the time 
step number /. Here, the 'consumption' function cf)'' and the function r' 
are defined recursively on the domains Df := [/ — 1,/] x M for / > 1 and 
depend on the solution v^''~^ and r'~^ defined at the previous time step. 
The functions 1 < i < n are defined by sums of source functions which 
are defined on different scales of magnitude. This ensures that the control 
of the local solutions v^'^''"'^^ does not interfer significantly with the control 
of r'''^^ . At the first time step I = 1 the dependence is restricted to the initial 
data set h. c) We ensure that the control function r and the function are 
globally Holder continuous with respect to space and time. This implies that 
classical arguments lead to classical C^'^-regularity of the velocity function 
V, and hence of the pressure. 



h G if' (M, TM) . 



(114) 



We prove 
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Theorem 3.1. Given any dimension n and a viscosity constant u > let 
h G C°° (M,TM) (or, equivalently satisfy 1^114^ for any s € Mj. Then there 
is a global classical solution 

V G C^'2([0,oo) X M,rM) (115) 

to the Navier-Stokes equation system in its Leray projection form. This 
solution is given in the form as a limit of a global scheme for v*" — r where 
is the globally Holder contiunous limit of the scheme above, and r is 
a bounded globally Holder continuous control function which is constructed 
within the scheme. 

Proof. We do the proof in four steps. 

1) In a first step we prove local convergence in a classical function space 
of the local scheme where r' = 0. This is the situation of the first 
time step Z = 1 of our scheme. Especially we determine a time step 
size. In a first substep we consider convergence of the spatially local 
iterations scheme with iteration index p. Note that in this substep 
the boundary conditions are independent of the iteration index. In 
a second substep we consider convergence of spatially global scheme 
with iteration index m. In this second substep we have a dependence 
of the boundary conditions on the iteration index m and we show that 
we have a time-local fixed point in a classical function space. 

2) In a second step we prove convergence of the time-local controlled 
scheme for v''''''' if the data r'~^(/ — 1, .) and v^'f'''~^(l — 1, .) are in 
C'^{M, TM). The scheme is similar to the local scheme of step 1) , but 
has a source term. Again, in a first substep we consider convergence 
of the spatially local iterations scheme with iteration index p at time 
step L In a second substep we consider convergence of spatially global 
scheme with iteration index m at time step I. In both substeps we 
assume regular bounded data. 

3) In a third step we show that the source term of step 2) is defined via 
a control function in our global scheme such that the growth of the 
absolute values of the components r\, 1 < i < n, and v^'^'\ 1 < i < n 
and the growth of the first order spatial derivatives of these functions 
is controlled independently of the time step number I > 1. 

4) In a fourth step we show the existence of a globally bounded classical 
solution G C^'^ ([0, oo) x M,TM) of the time-transformed incom- 
pressible Navier Stokes equation exists. It follows immediately that the 
global solution v G C^'^ ([0, oo) x M, TM) in original time coordinates 
exists. 
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3.1 step 1: proof of local existence of solutions at each time 
step (in absence of a control function) 

We emphasize that the time step size p is chosen generically, i.e. in this 
first step of the proof we shah choose p > such that local convergence is 
obtained in the scheme without control function. Similarly the bound C > 
of for the (modulus of) controlled velocities Vi + ri, I < i < n, the (modulus 
of ) control itself, and (the modulus of) their first derivatives will be chosen 
generically. We set I = 1 and r^'^^ = for all 1 < i < n and j £ J. The 
scheme for v'^'P'^'^^ = v^'^'^''^^ = v^'^''^^ is a scheme with iteration indices m 
and p. First for fixed m > 1 we consider for each p>l the equation 

dr t' l^q,k=l qk > dx^dx^ P l^k=l"k > dx^ 



+P22k=iv% 



(116) 

where for p = 1 we have -yP.'-'^J'^Cp-i) = yP,i,kjm{o) ._ yp,i,kj{m-i) ^ ^^^^ ^^^^^ 

the coupling term is independent of the iteration index p. For m = 1 we set 

^P,i,kj{m-i) ^ yp,i,kj{0) ._ yp,i-i,kj(^i _ gi^(,g / = 1 the latter function 

equals the initial data h^^ on Uj, where 1 < k < n denotes the component 
of h in local coordinates. However, keeping the notation vf'''~^'^^{l — 1,.) 
reminds that this is a part of a recursive scheme where an additional control 
function r' will be added later. At each iteration step p we defined 

yP,l,ijmp(^l - 1, .) = VP'^-^''^{1 - 1, .), (117) 

and for each j £ J we added a boundary condition 

^^''•^^■'"^l[z-i,z]xai/,(r,x) = vP'''^^>'^"'-'\t,x). (118) 

Note the double superscript kk which indicates the use of a partition of 
unity. This ensures that differentiability of the restriction of u^'' to 
[/ — l,i] X dUj in (jll8p . Hence we have no dependence on the boundary 
condition if we consider iteration with respect to the iteration index p. The 
solution v^'^'^^"^ (limit p t oo) may be represented in the form 

oo 

yP,l,ijm _ yp,ijml _|_ ^ ^yp,l,ijmp (119) 
p=2 

along with = vP'^''^"^P-vP'^''^"'^p^^\ For p = 1 we denote (^v'''''*^™^ = 

yP,i,ijmi _ yp,hijm{0) ^ yp,hijmi _ yP,i-i,ijm _ r^^As has the advantage of zero 
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initial conditions and zero boundary conditions for (5i;/''''*^™P for p>2. Fur- 
thermore the equation for is a linear parabolic equation with differ- 
entiable coefficients. The equation for Sv^'^''^^^'^ becomes 

+P Y.k-1 5v''{!''^'^^vP^^'^^'^''P~^^ 

= -P ELl vf'^'^^'^-^hvP^^^^i'^^P-l) (120) 

^ mt,i ^ ' 

— pS^. ^ (^■yr,l,ijm{p-2) ^j;;;-/^r,l,ijm{p-2)^ ^ 

with zero initial and boundary conditions which do not depend on the itera- 
tion index p. Note that the coupling term /oS'coup (v^'''*-'^™"-^), Vv'''''*-^("^~^)) 
disappears since it does not depend on the iteration index p. Hence the so- 
lution for the functions Sv'''''''^^'^^ for p > 1 has the representation 

fo'''^'^^--P(r,x) = I^^J^ 6s^^^^^^^^'\s,y)p'''^"'^P-'\T,x;s,y)dy^ (121) 
where p'>*-?™-(P"^) is the fundamental solution of 

dr ~ P l^q,k=l'^qk V^J dxqdxk ~P2^k=l"k \^ ) dx^ 

(122) 

+pZl^,vf''"'vP'^'>'^^(P-^^ = 0, 

and where '^•^j'^^^'^'^j^p ^\^iy) is an abreviation for the right side of (jl20p . 
It follows that for p > 1 we have 

^p,l,ijml^ ^p,l,ijmp ^ ^1,2 ([^ _ 1^/] X I^) . (123) 

In order to estimate (5f'''''*-'"*^(r, x) we may recall from classical theory that 
the fundamental solution in the Levy expansion on an arbirtrary 

domain — 1, Z] x Q form is given by 



p'''^^p{t, x; s, y) := iVi(r, x;s,y)+ j / N'^^t, x; a, OH^, s, y)dad^, 

(124) 



where for (a*-''-^'^^) defining the inverse of (ofj-^^ ) we have 



/det [a^^-.W] f a*.,W(^^_y^)(^. _y.) 

NA{r,x;s,y) = =n exp 



E • ,=1 4vrp(r - s) V MHr-s) 

(125) 
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and (p is a recursively defined function wliich is Holder continuous in x, i.e., 

oo 

4>{r, x; s,y) = J2 {LiN^A)m{r, x; s, y), (126) 

m=l 

along with the recursion 

{LiN^^)i{T,x;s,y) = LiN\{t,x; s,y) 



P 2^k=l "k ) 



(127) 



{LN'^U+i{t,x) := jlj^ {LN^^{T,x;a,0)^LN^Ai<',C;s,y)dadC. 

Note that in locally flat coordinates the entries of (a^-^'-^'^j) are locally con- 
stant on Uj such that for small p > the Levy expsnion is a kind of per- 
turbation of the leading term A'^^. Especially for second partial derivatives 
with respect to the variables x^ and Xq we have 

^^P,^jmp^^^ ^) ^ jr^^ J^^ Knt,COUp y)PkT^^~^^^^^ ^' 

+ II /r" ^A,k,qi'^^ CT, 0(/)(cr, ^; s, yUdad^dyds 



(128) 

and we observe that the integral involving the leading term involving N^j^ ^ ^, 
i.e. the second partial derivative of with respect to xj^ and Xqj is indeed 
a convolution. Hence we shift a derivative in the leading term writing 

^^P,lm^^^ a;) = /;^^ J^^ '^^inCcoup y^P'kT^'^^ ^' ^' y^'^^^^ 



= I^ik's^;^;^is,y){Nl,^^ir,x,s,y) 
+ II /m" ^A,k,q('^^ ^5 ^' O0(cr> ^; s> y))dad(,dyds 



(129) 

We can then estimate this directly or shift dervatives in the lower order 
terms inductively wich leads to a representation in terms of the adjoint of 
the fundamental solution. In any case inspecting the soorce term of the 
equation for (^u'^'^*"'™^, of lower derivatives, and observing that the time 
derivative can be locally estimated in terms of the spatial derivatives, for 
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small p > we observe that 

|5^P,Mjml|^ 2 < ^, (130) 

and 

2 < ^|5r;^'''^^''"(P-^)|i,2, (131) 
such that for this time-step size p > we have 

p=2 ^'^ 

(strictly less indeed). Hence for this time-step size p > the sequence 

^^r,p,l,jmp^^^^ (133) 

converge classically to a fixed point limit 

^r,p,l,ijm ._ ^^^^ ^r,p,l,ijmp ^ ^1,2 qq^ ^ jj ^ (-^34^ 

for all 1 < i < n and j G J, and in Banach space C^'^ ([/ — 1, x UjY For 
each m > 1 this fixed point limit solves the initial-boundary value prob- 
lem for stated in (fTHB . (fTHTl) . and (fT38]) below. Since we know the 
functions ti^'''*-?™- as the fixed point limits of subiteration steps with itera- 
tion index p for all 1 < i < n and all j G J, we an represent the functions 
yr,p,i,ijm terms of fundamental solutions pP^^''^^^ of the (nota bene nonlin- 
ear) equation 

dpP,i,ij"^ _ f,U j N| i9^pP''-'J'" _ v^n if,U I -.OpP^lAjj^ 

Or Pl^q,k=l^qk dxqdxk P 2^k=l"k dx^ 

(135) 

In order to write the known function u^'''*.'"^ in terms of this fundamental 
solution pP^^'"^!^ we note that for the functional sequence iy^'^'^'^^^) j^&i 
get another spatially global iteration scheme of local inital-boundary value 
problems 

dr P2^q,k=l"-qk ) dxqdxk P Ljk=l"k dxk 

+pY1=iv''^{!''''"vP'''^^^ (136) 

= ^^!nt'J(^'^^) + '^^co(lp"'(^'V^)' 

where for m = 1 we set = ^P^i^^m ■- yP,i-^Mj(j^ _ At each 

iteration step m we defined 

^P,i,ijm^l - 1, .) = vP^^-^^'\l - 1, .), (137) 
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and for each J € J we added a boundary condition (for {t,x) € [l — l,l]xdUj) 
v'^'^'nii-ifi.au.ir^x) = vP''''^'^^^-^\t,x). (138) 

For all (r, x) S Uj we have the representation 

vP^^''^'^{t,x) = f^^ vP'^~^''^{l - l,y)pP^^''^"'{T,x;0,y)dy 

+ III, s^^^^is, y)pP'^'^^^ir, x; s, y)dsdy (139) 

+ !i-i laUj '^bd(^' y)p'''''*^™(r, x; s, y)dSyds 
where we used the abbreviation 

mt.coup mt,i ^ ' ^ 

(140) 

P'^c^up (v^'''^-'('"-i), Vv^'''*-'(™-i)) , 

and where dSy denotes a surface element on dUj. The boundary relation 
reduces to an integral equation 

J^^ vP''-^''^{l - 1, y)pP''''i^{T, x- 0, y)dy 

+ /r-i h <nrSup(^' s, y)dsdy 

(141) 

+ !du, <^bd(*' y)p^'^''^'^{r, x- s, y)dSyds 

= EkeJ,y''''''^"'-'Hr,x) 

for the function 0^^. Again we may solve for n scalar initial boundary value 
problems for g C^'^ ([/ — 1,Z] x C/j) first, show that for small p > 

we have a small difference ti^'''*-?^ — ti^''~^'*-' (/ — 1, .) and then show that 
§yr,i,tjm _ yr,i,ijm _ ^r,i,tj{m-i) g^tigfies a Contraction 

\Sv'''^'^'^\i,2 < ^|5u''''^^'(''""^)|i,2. (142) 

This is done using the classical representations of initial-boundary value 
problems in terms of fundamental solutions as above in (jl39p and (|14ip . 
Note that we use the term for small p > in a generic sense here, i.e., we 
first determine a p such that we get the desired contraction for Sv^'^'^^"^^ 
with respect to the subiteration index p, and then we use this p in order 
to get another p which is smaller or equal such that the ()142p is satisfied. 
We do this here for one time step I, and the choice of the control function 
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willendure that it can be done independently of the time step number /. For 
a step size p > which is smah enough the sequences 

{-''''^nmen (143) 

converge to a classical limit 

for all 1 < i < n and j S J. Finally we set 

for all 1 < i < n and j S J. We choose C > such that 

|^r,p,l,*J|^ 2 < C. (146) 



3.2 step 2: Convergence of the spatially global and time- 
local controlled scheme 

We take the time-step size p > from the previous step in the generic 
sense, i.e., we may adjust this time-step size to a (possibly smaller) time- 
step size denoted also by p > such that we have time-local convergence 
in the presence of a local control function r'. Recursively, if v^'f'^~^''^^ , and 



l-l,ij 



are known, then we construct 



If, 1,13 _|_ M",i,^3 as described in 



and (|9Up above. Then we determine r''*-' via the initial-boundary value 
prolems described in ()82p with the initial and boundary conditions described 
there. Note that we determine the time-local function r''*-' via an iteration 
(although the equation is linear). The convergence of the iterative approxi- 
mations r''*-^™ ensures the existence of r' : — 1, /] xTM, i.e., the convergence 
of the respective boundary conditions for the functions r^'*-''" as m t oo. This 
involves all indices 1 < i < n and all j S J in order to make the scheme 
for the control function spatially global, of course. We shall choose C ^ 1 
(recall that the choice of C is generic) such that 



< C, and \v 



r,p,l—l,ij I 



and for all 1 < m < n 







Qyr,p,l-l,ij 


< C, and 







< c. 



Then we define 
and 



(147) 

(148) 

(149) 
(150) 
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Note that we put (f90)^150p on a different scale, i.e., we use a denominator 
in tlie definition of , sucli tliat tlie control of growth of the value 
function y'"'^'''*-? does not interfer with the control of growth of the control 
function r''*-' essentially. Then we compute the control function r''*-' as a limit 
of coupled initial boundary value problems, i.e., as a limit of the scheme 

dr P2^q,k=l"-qk \^)dxqdxk P2^k=l"k \^ ) dxk 

(151) 

-P^''^'"" - /'^ciu^'^^'""'^ Vv) , 

where 

- 1, .) = r^~^^i0m^i (;^52) 

and 

for (r, x) S [Z — 1, /] X Note the minus signs in front of pS'^^^^ ^'ii"^ i) g^y^^ 

r l—l J (m—l) 

in front of pS^^Q^^p . We choose them, and we choose them having in 

mind that we plug in the equation for r''*-' = lim^-i-oo r''*-'™ into the equation 
for ti^'/''''*-? for each 1 < i < n and j S J. For small p > we have convergence 

r''*J' = lim r''^J" G C^'^ ([Z - 1, /] x F") , (154) 

mtcxD 

and such that a global function r' : — 1,/] x Af — > TM is well-defined. 
For the spatially local iteration scheme for the sequence (?^^'^'''*-''"^)pgpj we 
defined the iteration scheme above pU4|) . We reorder the terms a little in 
order to prove contraction and write 
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Qyr,p,l,ijmp f^JJ, . g2 ^r,p,l ,ijmp „,„ l^f,Uf \ gyr,p,l ,ijmp 

Pl^j,k = l^jk dXjdXk Pl^k = l"k dTk 

= A},i3 _ pV" ^^^^iyr,P,l-l,kj _^ „ y " yr,p,l,kjm{p-l) 

-P ELi ^; J'^-i'^V'-i''^^- + p ELi vJ''^''"'^^-'\^'>'^ (155) 

-pS^r (v, Vv) + pS'^^^^^^ (v^sP,M.(— 1, Vv^.P.M.{— D) , 
with initial data 

^r,p,ZOT^Q^ ^ yr,p,l-l,ij _ (;^5g) 

and boundary data 

v'-^''^'''^"'P\[i.,,i]^9U,ir,x) = J2 v'^p^'^'^^^^-^\t,x). (157) 

In case m = 1 we define 

^^^'''''''^■("^"'^^I[^-M]x9C/,(t,x) = ^;^'^''-1'*'=(t,x) (158) 

for all (r, x) € [/ — 1, /] x dUj . We may write down representations of solutions 
for g^j^^ ^^r,p,l,ij[m-i)p terms of fundamental solutions of 

linear equations similar as in the first step of this proof. For a generic step 
size p > which is small enough the sequences 



^r,p,l,ijmp\ /-^5gN 



converge classically to a limit 



P'^''^"" := lim v'''f'^''^"'P G ([0, 1] x Uj) (160) 
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for all 1 < i < n and j € J. Finally, for the functional sequence {v^'''''''^'^"^) ^^j^ 
we get another spatially global iteration scheme 

dT Pl^j,k = l'^jk dXjdXk P l^k = l"k dXk 

En r,p,l—l,ij ;_i l-^ , v^n r,p,l,ij(m—l) i ki 
_|_ ^Y^n ,^l,ij^l,kj 



Pl^k=l^,k ^ + P l^k=l 



-P^lit'^' + /'^int, (v^'^'''^^(— 1), Vv^'^O 

-P^coupJ (v, Vv) +p5c^^up (v-.aMj(— i)^vv^''''''^^(— 1)) , 
where we rearranged the terms again, with initial data 

^r,p,l,ijm^^ - 1, .) = V'''P'^-^''^{1 - 1, .), (162) 

and with boundary data 



r,p,l,ijm 



V—\[l-i,l]xdU,{r,x) = ^^''''^•^^'=(™-i)(r,x). (163) 

Again for m = 1 we may use time-constant data 

v'''''''%-i,l]xdU,ir,x) = v^-'P'l-^'^\l - 1, .). 
For small p > we get a classical limit 

^r,p,l,ij ._ ^-^^ ^r,p,Mim ^ (-1,2 Qg, 1] x Uj) (164) 

m.^oo 

for all 1 < i < n and all j S J, and such that a global function v^'^'' : 
[/ - 1, /] X M ^ TM is well-defined. 

3.3 step 3: Control of the growth of the functions r' and v^ '''' 
In this step we shall adjust the generic time-step size p > such that 

P^^y (165) 

Here the constant C > is also generic (greater or equal to C > chosen 
at the first time step and may be assume to be greater w.l.o.g.). The reason 
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for a choice of type ()165p is our definition of the source function cj}''^^ which 
determines the growth of the controlled Navier-Stokes equations for ti^'/''''*-? . 
We noted that we define 

4>''' - ^, (166) 

a behaviour which is determined by one summand of cj}''^^ which controls 
the growth of ti^'/''''*-? . We have proved that r''*-' and v'^'P'^''^^ exist for all 
1 < z < n, all j € J and all / > 1 as functions in C^'^ {[1 — 1,1) x Uj) provided 
r'^^'*'' (Z — 1, .) and v'^'P'^"^'^^ {1 — 1, .) exist as functions with continuous spatial 
derivatives of first order. Moreover, by imposing boundary conditions for the 
approximations r'"*-^"*, we have shown that the local solutions of the linear 
equations for the control function r''*-' = lim^-i-oo r match on common 
domains Uj n Uk for all j,k € J, whenever the latter intersection is not 
empty. Similar for v^'P''"''^^. This defines for all / > 1 spatially global and 
time-local control functions 

: [1-1,1] X M ^TM, (167) 

provided that r'~^ and v^'^''~^ exist and are sufficiently regular. Similarly, 
the analyis of the time-local iteration of our scheme in the first and second 
step of this proof has shown that there exists a spatially global function 

V^P^^ -.[1-1,1] X M ^TM, (168) 

whenever we have well-defined initial data r'~^(/ — 1, .) and v'^'P'''~^{1 — 1, .) 
which are in C^. Next we have to show that r''*-' and v'^'P'^'^^ are uniformwith 
respect to time). First we show that the bound for r''*-^ and its first order 
spatial derivatives are preserved for each time step I. Then we show that the 
bound for v^'P'^''^^ and its first order spatial derivatives are preserved as well 
for all time steps /, where we use the fact that the source function cf}^'"^'^^ is 
defined on a different scale of magnitude, let us show first that the growth 
of r''*-^ is bounded. More precisely, let us show that the absolute value of r''*-' 
and the modulus of first spatial derivatives of the functions r''*-' are bounded 
by a constant C > which is independent of the time step number / (and 
trivially independent of the index 1 <i <n and j € J- note that J is a set 
of finite cardinality) . First we consider a representation of a solution of the 
linear parabolic equation ()82p with initial data r''*-' (Z — 1, .) = r'~^'*-' (/ — 1, .) 
in terms of the fundamental solution F'^'''*-' of the equation 



ar'-'J ^Y^n nf^U r ^&^r^ _ sr"^ ^/.^/'^^ Qr-'-'J 
dr P2^q,k=l"'qk \-^)dxqdxk P2^k=l"k ) dxk 



(169) 



on the domain \l — l,l]xUj. Note that these function r''*-' where constructed 
by an iteration with iteration index m which ensures that they all match on 
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nonempty intersections Uj D Uk for j, k ^ J. Then a solution of the equation 
([82]) for r^'*-' with initial data r''*-' and boundary data r'''"^'*-' on [/ — 1, x dUj , 
i.e. 

r^M,ij^^^ ^) = x) = Y^ r^'^'^Hr, x), (170) 

has the representation 

r'.*i(r, x) = /j^^^ r'-i'^J' (/ - 1, y) r'''^^ (r, x-l-l, y)dy+ 

(171) 

- T-i /l/, P^hi''' (v, Vv) (s, y)r-.M.(^, z _ \^y)dyds 

- Vi-x lu, P^coup^" (v, Vv) (s, y)r'''*^(T, x; / - 1, y)dyds 
+Pli-Jdu, <t>'^'{s,y)r^^^''HT,x;l- l,y)dSyds 

=: r{T,x) + Jqij. (/'|,']^(s,y)r'''^''J(r,x;/ - l,y)dSyds, 

where in the latter term dSy is the surface element on dUj. The function 
^bd boundary condition and, hence, is determined via the integral 

equation 

r / '^l'r(^'?/)r''''''(^'^;^-l'y)rfVs = ^(r,x), (172) 
Jl-1 JdUj 

where 

R{t,x) := Y.k&J, r^''^^{r,x) - r{T,x). (173) 

Invoking the maximum principle we observe that the maximum value of 
on Uj is not exceeded by the first term of ()17ip -since classical deriva- 
tives exist apply a maximum principle iteratively (with respect to the iter- 
ation index m) to r^'*-^™ — r^^^; where 

^max := maXjGjrmax := max^gj maxr^ax - 1> a^)- (174) 

Similarly for — r''*-' . Furthermore, for small p > the source term is the 
only source term without a factor p and dominates all other source terms. 
Hence we get 

max |r'''^(r,x)| < maxr^ii^^^(/ - l,x) < C, (175) 

(r,a;)G[i-l,i]xC/j x<^Uj 
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inductively, i.e., if this bound for r^'*-' holds for time step / — 1 then it holds 
for time step / (for all 1 < i < n and all j € J. This p > can be chosen in 
the form ^ 

P= / N^3 . (176) 

where C > is the bound with respect to the |.|i^2-iiorm of the previous 
time step I — 1 and c(n) is a constant depending on dimension n and on the 
number of terms appearing on the right side of the equations for r''*-' and 
yi',p,i,ij respectively. Note that we consider here an iteration with respect 
to the iteration index m because we construct a spatially global control 
function r' at each time step I. Concerning the bounded first derivative of 
r''*-' we invoke the Levy expansion, locally flat coordinates, and small time 
step size /) > 0. 
Hence we get 

max |^r''*^(r,rE)| < C, (177) 
{T,x)(^[l-l,l]y.Uj ' OXk 

inductively for all 1 < A; < n, i.e., the estimate (|187p holds if it holds for 
ri-^,i3(i _ 1^.). 

Next we look at the scheme for We have proved that time-local 

limit t;'''^'^'*-' = limmtoo.ptoo v^'^''''^^"^^ exists and satisfies for all 1 < i < n and 
all j G J the equation 

dr Pl^q,k=l^qk \^ > dxgdxk P2^k=l"k \^ ) dx^ 

- dr P 2^q,k=l°'qk dxqdxk P 2^k=l'^k l^i&fc \^'^) 

+P ELi rf + p ELi vf''^r'^^^ 

+P ELi rffr^^^ + pS^^^ (v, Vv) + pf^l;^^ (v, Vv) , 
on the domain [/ — 1] x Uj with initial data 

yr,p,l,ij^l _!,.)= i;-'P''-l'^i(Z - 1, .), (179) 

and with the boundary condition 

^^''"''''^l[o,i]xai/,(r,x) = v^'''''''^\t,x). (180) 
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Next we plug in the limit of ()82p . i.e., 

dr P Z-^q,k=l"'qk \-^)dxqdxk Pl-,k=l"k y-^> dx^ 

= -pEk=i rfr^-^''^^ -pEI=i rjf t;^'^''-!'^-^- - pELi vj^'^'^'^r^-^'''^ 
~'^'^int~i ~ P^cm/' (v, Vv) , 



into the right side of p78p . This leads to the equation 

dr P 2^q,k=l^qk \-^) dxgdx^ P l^k=l"k dx^ 



(181) 



P 



-^^[nt"''' + ^^[nt'. - /^^coup^^ (V, Vv) + pS[i^ (v, Vv) . 

(182) 

Note that for small p > the source term (p^'^^ is the dominant term on the 
right side of the latter equation. We may assume that C > is large enough 
such that the source 

<^^.^.*i(x) := -^A^-^)^iJ(x) (183) 

for r''*-' dominates the source 

</.^'''^i(x) := -^^;'^.P.('-i).^i(x) (184) 

for ti^'^'''*.? whenever r'"^'*-' is large. Nevertheless we may choose p > small 
enough such that the source (|184p is small enough in order to control the 
growth of the function v^''''^'^^ (and its first spatial derivatives), whenever 
the control function r^~^'^^{l — 1,.) is small. Since v^'^'^'^^ is known after 
time step I we have a representation of the initial-boundary value problems 
for v^'P'^'^^ in terms of the fundamental solution F"''''^-^ of the equation 

q.k=l « k=l k=l 

(185) 

For C > large enough the additional source terms (jf^^'"^^ are small compared 
to the source term such that a bound for the control of growth of 

the control function r is determined by (j/'^^'^K For the growth of u'"''''''*-' 
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is controlled by the source function (ff'^''^^ in all regions where the modulus 
of r''*-^ or the modulus of its first derivatives exceed a certain level. The 
small step-size p ensures, that the diffusion effect is small and all terms with 
factor p are small, such that the source term dominates the growth behavior. 
Otherwise, in regions where the control functions r^'^^ are small the smaller 
source (/("'^'^J still dominates the growth behavior of the controlled velocities 
yr,p,i,ij g[T^QQ ga^]2 choose the time-srep size p > small compared to the 
factor A similar argument as above in the case of r''*-' then leads to the 
conclusion that 

for all 1 < 2 < n and all j G J and all / > 1, and 

max \-^v''P'^''^{t,x)\ <C, (187) 

{T,x)G[l-l,l]xUj OXk ' 

for all 1 < i, A: < n and all j J and all / > 1, inductively. 



3.4 step 4: Global existence of classical solutions and v 

Now we have proved that for all 1 < i < n, all j € J, and for all / > 1 we 
have 

^c^^[l-l,l]x Uj) , and |r''*^|5 < C (188) 

and 

^r,p,i,ij ^ ^5 _ 1^ /] X Uj) and \v'-'P'^''^\s < C (189) 

for a constant C > which is independent of the time-step number /. Indeed 
we have more regularity with respect to the spatial variables and even with 
respect to the time- variable r (transformed time) we have classical differen- 
tiability except at the points r = 1, 2, • • • , i.e., where r is a natural number. 
The non-differentiability in a classical sense with respect to time at these 
points is due to the fact that the source functions of the equations for 
r''*-' are locally constant with respect to the time variable over time [Z — 1, /) 
and we have bounded jumps from to 0''*-' at time r = / — 1 in general. 

However these source terms appear as a time integral in the representation 
for r''*-^ and for v^'f^'^'^^ (or its first approximation), and this leads to the con- 
clusion that r''*-^ and v'^'^'''''^^ are Holder continuous across the time points 
r = I for all time step numbers Z > 1. Now we immediately get for all 
1 < i < n, all j G J, and all / > 1 that 

^P,i,ij g fjS _ 1^ /] X Uj) and \vP^^''^\s < 2C, (190) 
with the same C > independent of / > 1, and where 
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Hence 

yP,ij g ([0, oo) X Uj) and \vP''^\s < 2C, (192) 

where we recall that f^'*-' is the function which equals ti^'''*-? if restricted to 
[1 — 1,1) X Uj for alH > 1. The next observation from the argument of the 
preceding steps is that the first spatial derivatives of r''*-' and of ti^'/''''*-' exist 
continuously for alH > 1 for all 1 < f < n and all j G J. Therefore, we have 
for all I > 1, for all 1 < i < n, all j £ J, and all 1 < A; < n 

v'^^'^ G C([0,oo) X Uj) and |';^''''-''^'lco(I7-) < 2C- (193) 

This holds for all j € J. Then looking at ()3ip we observe that the first order 
coefficient v^^ satisfies 

v^^ eC^ {[0,oo) xUj), (194) 

i.e. the first order coefficients are Holder continous, and this holds also for 
the intergal terms on the right side (invoking regularity results for Poisson 
equations). Hence, from classical theory of linear parabolic equations we get 
for all 1 < i < n and all j € J 

yP''^ £C^'^{[0,oo)xUj). (195) 

It follows immendiately that this holds also without time dilatation, i.e., for 
all 1 < i < n and all j € J 

v'^ €C^^^ {[0,oo) xUj) , (196) 

where we recall that v^^ {t, .) = vP'^^{t, .) along with t = pr. Hence, for the 
restriction of the velocity components f * to Uj we have 

v%^eC^'^{[0,oo)xUj) (197) 

for all 1 < i < n and all j G J . It follows that 

V G C^'2([0,oo) X M,rM) , (198) 

and our argument is finished. □ 

Finally we note that the existence of a classical solution implies that the 
solution is smooth, i.e., we have 

Corollary 3.2. For s smooth Riemannian manifold M a viscosity constant 
u > and data h G C°°{M,TM) we have 

V G ([0, oo) X M, TM) (199) 
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Proof. We know v G C^'^ ([0, oo) X M, TM) from our main theorem. This 
means that locahy with v^^ known on Uj for ah j £ J we may look at the 
local incompressible Navier-Stokes equation as an uniform parabolic equa- 
tion which is linear in terms of v^^ and its first spatial derivatives, i.e. we 
look at ()3ip in the form 

-w l^j,k=i"'jk y^'dxjdxk \l^k=i'^k y^) ^ )'d^ 

(200) 

= /, 

where 

f = 4t, + ^coup,i ' (201) 

and we observe form the definition of /, our main theorem, and standard 
Schauder estimates that the first order spatial derivatives of the right side f^i 
are Holder continuous. Then we may iteratively apply the classical theorem 

Theorem 3.3. Assume that for all multiindices a with \a\ < m we know 
that 

D^i,, D^k, D^c (202) 

are Holder continuous in a domain D = (0, oo) xQ for some domain Q, C M". 
If u is a solution of 

in D, then 

D^u, DtD^u (204) 

exist for < |a| <m + 2 and < <m and are all Holder continuous in 
D. 

This leads to f*-' G ([0, oo) x Uj) for all j G J and then to to global 

smoothness as stated. □ 

The explanation give here for a global scheme of the incompressible 
Navier-Stokes equation uses strict and uniform parabolicity. It is clear that 
this assumption cannot be removed for we know that solution of Euler- 
equations (with viscosity v = Q) may blow up in finite time. However, the 
method considered here may be extended to systems where the second order 
coefficient functions satisfy a Hormander condition. The reason is that in 
this case Gaussian density estimates of Stroock-Kusuoka type seem to be suf- 
ficient. Each time step then involves the solution of problems on [0, oo) x M 
of the form 

(205) 

w(0,x) = f{x). 
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for m vector fields which look locally like 

" d 

Vi = Y.v,,{x) — , (206) 

j=i J 

where < i < m. For equation on manifolds the Hormander condition has 
to be rephrased in local charts, of course. However, independence of the 
chart and well-definiteness is straightforward. Denoting the vecor fields Vi 
in a given chart on an open set 0, C by Vi again (same name) we may 
say that the Hormander condition is satisfied at x € if 

:= spanj Vi{x), [Vj, Vk] (x), 

(207) 

[[Vj,Vk],Vi]{x),---\l<i<m, 0<j,k,l---<m^ 

holds. We may say that the Hormander condition is satisfied on ^1 if ()207p 
is satisfied for all x G Note that second order equations of type ()205p 
correspond to diffusion processes X which have a local representation on 
domains in components, and satisfy stochastic ODEs. In the framework 
of Malliavin calculus it was proved that 

Theorem 3.4. Let the assumption of be satisfied and let T > 0. Then 
the law of the diffusion process X ( corresponding to the second order equation 
1205\) in the Feynman-Kac sense) exists on a domain C is absolutely 
continuous with respect to the Lebesgue measure, and the density p exists 
and is smooth, i.e., on a domain 0, C R" we have 

p : (0, T] X X 17 ^ R G C°° ((0, T] x Q x Q) . (208) 

Moreover, for each nonnegative natural number j, and multiindices a, f3 there 
are increasing functions of time 

A-,a,/3,^j,a,/3: [0,r] ^R, (209) 



and functions 
such that 



nj,a,p, nij^o^^fi : N X X N'^ ^ N, (210) 
P{t,x,y) < , exp { -Bj^aA^^) 



dp dx"' dyf^^ ' ' - i"j.">/3 ^ V -"""^^ ' t 

(211) 

Moreover, all functions \209\) and 121 0\) depend on the level of iteration of 
Lie-bracket iteration at which the Hormander condition becomes true. 

These density estimates fit in our scheme and may lead to generalisations. 
Note that polynomial growth factor (1 + x)"^^'"''' in (j21ip is no obstacle 
since we work on compact manifolds. However, there are some additional 
obstacles, and will be investigated further in a subsequent paper. 
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